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ABSTRACT 

A  comprehensive  methodology  is  presented  for  designing  instrumented  oceanographic  surface 
moorings  for  harsh  weather  environments.  In  the  pasty  the  design  of  oceanographic  mooHngs  was  based  on 
static  analysis  as  most  systems  were  deployed  m  regions  where  wave  forcing  was  small.  But  now,  many 
surface  moorings  are  being  placed  :n  ocean  environments  where  the  predominant  forcing  is  from  waves 
and  where  the  major  cause  of  failure  is  cyclic  fatigue.  HerCy  dynamic  analysis  becomes  as  important  as 
the  static  calculations.  In  this  paper,  we  present  the  equations  and  numerical  solution  techniques  for  the 
statics  and  dynamics  of  a  buoy/cable  system  with  attached  instruments.  Hydrodynamic  coefficients  and 
wave  exciting  forces  of  diffierent  shaped  oceanographic  buoys  are  given  along  with  the  material  parame¬ 
ters  for  the  components  and  instruments  that  make  up  the  mooring  line.  We  show  how  to  combine  the 
results  of  the  static  and  dynamic  analysis  with  data  from  laboratory  strength  and  fatigue  tests  of  mooring 
components  to  predict  if  the  mooring  will  survive  the  deployment.  We  present  an  analytical  model  of  the 
dynamics  of  the  mooring  system  that  gives  almost  identical  results  with  the  numemcal  simulation  and 
can  be  used  during  preliminary  design.  A  comprehensive  design  example  using  the  1993  Arabian  Sea 
surface  mooring  is  presented  to  illustrate  the  design  procedure. 


INTRODUCTION 

Oceanographers  rely  hea\dly  on  single-point 
surface  moorings  with  attached  instruments  to  make 
measurements  of  current,  temperature,  conductivity, 
light  transmissibilit}’*,  fluorescence,  and  other  scientific 
quantities  in  the  water-column.  No  other  technique 
exists  that  can  provide  these  measurements  close  to 
the  surface  and  simultaneous  with  meteorological  mea¬ 
surements  above  the  surface.  This  is  important  to  sci¬ 
entists  who  are  studying  the  relationship  between  at¬ 
mospheric  forcing  auid  ocean  response.  Surface  moor¬ 
ings  also  allow  scientists  to  receive  d^ta  in  real  time 
via  satellites. 

The  history'  of  oceanographic  surface  moor¬ 
ings  at  Woods  Hole  Oceanographic  Institution 
(WHOI)  dates  back  35  years.  At  that  time,  oceanog¬ 
raphers  envisioned  a  mooring  with  attached  current 
meters  that  would  be  low  cost  and  could  be  deployed 
easily  from  a  small  ship  in  water  depths  that  ranged 
from  2000  m  to  over  5000  m.  Their  initial  deploy¬ 


ments,  however,  showed  how  diincult  it  was  to  main¬ 
tain  a  small  surface  mooring  for  periods  greater  than 
rvo  months  (Richardson  et  al.j  1963).  Failure  rates 
were  high  with  the  principle  causes  being  kinks  form¬ 
ing  in  wire-rope  mooring  line,  corrosion  of  the  wire 
rope,  fish  biting  through  synthetic  mooring  line,  and 
breaking  of  steel  mooring  hardware. 

A  major  engineering  study  in  the  late  1960's 
and  early  1970 ^s  (Berteaux  and  Walden,  1970)  ad- 
*■  dressed  these  failures  and  led  to  the  design  of  the 
WHOI  taut  surface  mooring  (Fig.  la)  which  became  a 
protovyipe  for  deploying  instrumentation  in  deep  ocean 
environments  (Berteaux,  1976).  The  design  consists 
of  a  combination  of  chain  and  wire  rope  in  the  upper 
1500-2000  meters  to  protect  against  fish  bites  and  ny¬ 
lon  rope  in  the  lower  part  of  the  mooring  to  give  com- 
pliatnce.  The  wire  rope  uses  a  3x19  torque-balanced 
construction  to  prevent  kinking.  It  is  galvanized  and 
rubber  coated  to  prevent  corrosion.  The  ends  of  the 
wire  rope,  where  it  attaches  to  the  instrument  cages, 
are  swaged  into  stainless-steel  fittings  and  covered  with 
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Figure  1:  (a)  Taur  oceanographic  surface  mooring  con¬ 
sisting  of  an  upper  string  of  instruments,  '^ire  rope, 
and  a  lower  section  of  nylon  rope  that  allows  the  moor¬ 
ing  to  stretch  in  strong  currents,  (b)  Inverse-catenary 
mooring  that  has  an  extra  section  of  polypropylene 
rope  near  the  bottom  to  increase  scope  and  reduce 
mean  tension  in  currents. 

water-tight  polyurethane  boots.  At  the  bottom  of  the 
mooring  is  a  recovery  system  that  uses  an  acousti¬ 
cally  controlled  release  placed  under  a  string  of  pos¬ 
itively  buoyant  17-inch  diameter  glass  spheres.  When 
released,  the  glass  spheres  float  to  the  surface  allowing 
for  retrie\"al  of  all  the  instruments.  The  anchor  and 
anchor  tackle  remain  on  the  bottom.  The  size  of  the 
anchor  is  chosen  so  that  it  will  lift  off  the  bottom  (al¬ 
lowing  the  mooring  to  drift)  before  the  surface  buoy  is 
pulled  underwater.  The  advantages  of  the  taut  moor¬ 
ing  are  that  it  has  a  smadl  watch  drcle  and  the  attached 
instruments  are  kept  nearly  vertical.  The  disadvantage  .. 
is  that,  in  strong  currents,  the  mean  tensions  can  be¬ 
come  large  enough  so  that  the  mooring  will  drift.  Also, 
as  we  show  later,  the  taut  mooring  cam  experience  high 
dTOamic  in  large  sea  states. 

The  instruments  are  attached  in-line  on  the 
mooring,  usually  being  concentrated  in  the  upper 
chain/ wire-rope  section  of  the  mooring  line.  On  some 
moorings,  there  have  been  as  many  as  15  scientific 
and  engineering  instruments.  On  WHOI  moorings,  the 
workhorse  instrument  is  the  Vector  Measuring  Current 


Figure  2:  Instrumentation  for  oceanographic  sur¬ 
face  moorings,  (a)  Vector  Measuring  Current  Meter 
(WICM)  for  measuring  current  speed  and  direction, 
(b)  SEACAT  for  measuring  conductivity  and  temper¬ 
ature,  (c)  Acoustic  Doppler  Current  Profiler  (.ADCP) 
for  measuring  current  profile. 


Meter  (VMCM)  (Weller  and  Davis,  1980)  and  more, 
recently,  an  amoustic  doppler  current  profiler  (ADCP). 
In  both  of  these,  the  electronics  are  placed  in  pressure- 
resistant  housings  and  suspended  in  a  cage  made  of 
3/4-inch  diameter  stainless-steel  rod.  Other  t3/pes  of 
instruments  have  been  placed  in  cylindrical  pressure 
housings  and  attached  to  a  single  1-1/4-inch  diameter 
stauniess-steel  rod  which  acts  as  the  strength  member 
(Fig.  2).  The  steel  cages  amd  rods  are  attached  to  the 
chain  or  the  swage  fittings  at  the  end  of  the  wire  rope 
using  galvanized  steel  shackles  and  links  (Fig.  3). 

A  variation  of  the  taut  mooring  is  the  inverse- 
catenary  mooring  (Fig.  lb).  The  design  was  orig¬ 
inally  proposed  in  the  late  1950's  for  mooring  the 
Navy  Oceanographic  Meteorological  Automatic  De¬ 
vice  (NOMAD)  buoy  (Berteaux,  1976).  Since  then, 
it  has  been  used  extensively  by  the  National  Data 
Buoy  Center  (Canada  and  May,  1985)  for  their  offshore 
weather  moorings.  More  recently,  the  inverse-catenary 
mooring  has  been  adapted  for  deploying  oceanographic 
instrumentation  in  severe  current  and  wave  environ¬ 
ments.  In  this  design,  a  piece  of  positively  buoyant 
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Figure  3:  Instruments  attached  in-line  on  the  mooring 
with  shackles  and  oval  end  links. 


polypropylene  rope  is  inserted  near  the  bottomjust 
above  the  glass-bail  recovery  systemrin  order  td  in¬ 
crease  scope  and  at  the  same  time  prevent  tangling 
of  excess  line  around  the  glass-ball  recovery  sys¬ 
tem.  The  extra  scope  helps  to  reduce  static  ten¬ 
sions  through  changes  in  the  geometry  of  the  mooring 
line  rather  than  through  stretching.  However,  hydro- 
dynamic  damping  constrains  geometrical  changes  at 
higher  wave  frequencies  (Triantafyllou,  1987)  so  that 


Figure  4;  Fatigue  failure  of  5/8-inch  pear-shaped  sling 
link  on  the  1989  MLML  mooring. 


the  d\mamic  response  is  similar  to  that  of  the  taut 
oceanographic  surface  mooring. 

The  success  of  WHOI  surface  moorings  dur¬ 
ing  oceanographic  experiments  of  the  IQSO's  can  be 
attributed  to  a  systematic  design  procedure  that  was 
developed  for  analyzing  the  static  configuration  of  a 
mooring  in  currents  (Berteaux,  1991).  The  success 
can  also  be  attributed  to  the  fact  that  these  early  de¬ 
ployments  (Briscoe  and  Weller,  1984;  Weller  et  a/., 
1990)  took  place  in  environments  that,  except  for  iso¬ 
lated  storms,  experienced  small  sea  states.  However, 
there  has  been  a  push  in  the  past  few  years  to  deploy 
oceanographic  surface  moorings  in  harsh- we  at  her  en¬ 
vironments.  Areas  in  the  North  .\tlantic  and  in  the 
Arabian  Sea  which  experience  high  winds  and  large 
seas  a  significant  portion  of  the  tune  have  been  sites 
of  recent  moorings.  One  mooring  in  particular,  the 
1989  Marine- Light j  Mixed-Layer  taut  surface 

mooring  deployed  400  km  south  of  Iceland  (Pleud- 
demann  et  of.,  1995),  highlighted  the  need  for  in¬ 
corporating  d}Tiamic  analysis  into  the  design  process. 
This  mooring  parted  after  70  days  when  a  5/8-inch 
galvanized-steel  sling  link  placed  at  the  top  of  the 
mooring  line  failed  in  fatigue  (Fig.  4).  This  is  in  spite 
of  the  static  analysis  that  predicted  the  component 
could  easily  withstand  the  maximum  expected  loads. 

What  was  missing  in  the  analysis  of  the  1989 
MLML  mooring  was  a  procedure  for  predicting  dy¬ 
namic  tensions  and  fatigue  in  mooring  line  compo¬ 
nents.  Because  of  the  compliance  built  into  taut  and 
inverse-cat enarv*  moorings  through  the  s>'nthetic  rope, 
the  instrument  string  is  free  to  accelerate  up  and  down 
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Design  Flow  Chart  for 
Oceanographic  Moorings 


Figure  5:  Flow  chan  showing  how  stadc  and  d}Tiamic 
analysis  are  combined  into  ultimare-strengch  and  fa¬ 
tigue  analysis. 


nearly  in  phase  wii:h  the  buoy.  Depending  on  the  to¬ 
tal  mass  of  the  instruments  and  the  amplitude  and 
frequency  of  the  waves,  large  dynamic  tensions  can  re- 
siilt  at  the  top  of  the  mooring  line.  These  oscillatory 
loads  can  be  of  similar  magnitude  to  the  static  loading, 
and  they  are  applied  over  and  over  again  as  many  as 
10  million  times  during  a  one-year  deployment.  Even 
in  relatively  benign  wave  environments  where  the  goal^ 
is  to  madntain  a  mooring  for  up  to  one  year,  dynamic 
tensions  can  play  a  major  role  in  the  survival.  This  was 
made  clear  when  two  moorings  of  the  1993  Sub4uction 
experiment,  deployed  1200  km  west^of  the  Canary  Is¬ 
lands,  parted  when  titanium  current  meter  cages  failed 
in  fatigue  at  weld  joints. 

The  goal  of  this  paper  is  to  present  a  com¬ 
prehensive  methodology  for  designing  instrumented 
oceanographic  surface  moorings  for  harsh  weather  en¬ 
vironments  that  takes  into  account  the  effects  of  both 
static  and  dynamic  loading.  The  important  point  here 


is  that  by  incorporating  dynamic  loading  we  are  able 
to  estimate  fatigue  life  and  to  make  more  accurate 
predictions  of  extreme  tensions  which  are  a  combi¬ 
nation  of  static  and  dynamic  tension.  Fig.  5  gives 
a  flow  diagram  of  how  we  view  the  ::.cign  process. 
The  main  steps  are  characterizing  the  .oean  environ¬ 
ment  in  terms  of  wnnd,  waves,  and  current.  This  is 
followed  by  static  and  dynamic  analysis  for  predicting 
the  different  loads  that  a  mooring  component  is  likely 
to  experience.  Finally,  using  the  material  properties 
along  with  safety-  factors,  appropriate  components  are 
selected  that  can  resist  both  the  extreme  loads  which 
cause  ultimate  failures  and  the  cyclic  loads  which  are 
responsible  for  fatigue  failures. 

While  the  offshore  industry  is  well  versed  in 
many  of  the  areas  that  we  address  in  this  paper,  the 
oceanographic  community  has  ignored  these  concepts 
in  the  design  of  their  moorings.  One  of  the  reason’s 
for  this  is  today's  oceanographic  engineers  are  required 
to  perform  a  wide  'v-ariety  of  jobs.  Besides  designing 
moorings,  they  design  and  build  instruments  that  go 
on  the  moorings,  and  they  go  to  sea  to  supervise  de¬ 
ployment  and  recovery.  Because  of  these  demeuids,  it  is 
very  difficuit  for  them  to  specialize  in  dynamic  analysis 
aind  fatigue  analysis.  And  because  of  costs,  it  is  even 
more  difficuit  to  hire  outside  consultants.  With  this 
paper,  we  present  the  theory  for  doing  sophisticated 
mooring  analy-sis.  the  hydrodynamic  atnd  material  con¬ 
stants  for  t3v-pical  buoy  and  mooring  components,  and 
examples  that  show  engineers  how  to  apply  the  con¬ 
cepts  to  their  designs. 

Beyond  this,  w-e  tailor  the  analysis  to  the  spe¬ 
cific  structure  of  an  instrumented  oceanographic  moor¬ 
ing.  Because  of  the  small  relative  size  of  the  buoy, 
interaction  betwreen  the  mooring  cable  and  the  buoy 
must  be  taken  into  account.  For  large  moored  offshore 
structures,  cable-structure  interactions  can  be  simpli¬ 
fied  by  using  a  quasi-static  approach  which  assumes 
that  the  dymamics  of  the  cable  do  not  effect  the  mo¬ 
tion  of  the  buoy-. 

In  the  next  section,  we  reriew  procedures  for 
calculating  the  static  configuration  of  the  mooring  and 
present  numerical  and  analytical  tools  for  predicting 
dynamic  loads.  This  includes  a  presentation  of  meth¬ 
ods  for  calculating  hydrodynamic  coefficients  of  small 
axisymmetric  buoy-s  which  are  needed  for  the  dynamic 
calculations.  Sample  calculations  aire  presented  that 
show  that  the  dynamic  tensions  in  the  mooring  line  are 
controlled  primarily  by  the  heave  motion  of  the  buoy 
and  the  longitudinal  motion  of  the  mooring.  We  use 
this  result  in  Section  3  to  construct  a  simple  maiyt- 
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icai  model  for  predicting  djmamic  tensions  in  surface 
moorings  that  gives  comparable  results  to  the  numer¬ 
ical  method  and  matches  full-scaie  experimental  data. 
In  Section  4,  we  show  how  to  combine  results  from 
static  and  d\mamic  analysis  to  predict  ultimate  fail¬ 
ures  and  fatigue  failures.  A  design  example  is  given  in 
Section  5  followed  by  conclusions. 

STATIC  AND  DYNAMIC  CALCULATIONS 

The  numerical  calculations  for  the  statics  and 
dynamics,  which  we  present  in  this  section,  are  based 
on  the  equations  of  motion  of  the  cable,  attached  in¬ 
struments,  and  the  surface  buoy.  We  assume  that  all 
the  motion  is  confined  to  a  single  plane  so  that  a  two- 
dimensional  formulation  applies.  For  the  dynamics, 
we  linearize  motion  about  the  nonlinear  static  config¬ 
uration  and  perform  frequency  domain  calculations. 
Our  mathematical  formulation  follows  the  procedures 
described  in  Goodman  et  aL  (1972).  The  numerical 
calculations  for  the  statics  are  identical  to  Goodman  et 
d,  (1972)  while  the  dynamics  are  modified  to  include 
nonlinear  drag  effects  through  an  equivalent  lineariza¬ 
tion  procedure  (Triantafyllou  et  a/.,  1987).  We  also 
present  a  novel  scheme  that  reduces  the  solution  to  a 
simple  closed  form  matrix  equation. 

Static  Calculations 

We  start  with  the  equations  for  the  static 
equilibrium  of  a  cable  element  of  unstretched  length 
ds,  weight  in  water  per  unit  unstretched  length  Wq, 
stretched  cross-sectional  area  .4,  and  Young's  modu¬ 
lus  E,  Using  the  derivation  and  nomenclature  given 
in  Triantafrdlou  et  d.  (1987)  we  can  wTite: 

=  Wo  sin  -  Fto(l -i- Co)  (1) 


dig  1 

where  To  is  the  effective  static  tension, (Goodman* and 
Breslin,  1976),  po  is  the  static  equilibrium  angle  of  the 
cable,  and  to  —  TojEA  is  the  static  strain.  The  co¬ 
ordinate  system  is  taken  along  the  cable  with  position 
3  and  unit  vectors  i  emd  n  representing  the  tangential 
(along  the  axis  of  the  cable)  and  normad  (perpendicular 
to  the  axis  of  the  cable)  directions  (Fig.  6).  The  fluid 
forces,  resolved  along  the  coordinate  axis,  are  given  in 


z 
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Figure  6:  Cable  coordinate  system. 

(Breslin,  1974).  The  fluid  force  tangential  to  the  cable 
axis  is; 

FtoiX^eo)  =  ^pCot^DoU cos  (po\U  cos  doWl  -r 

(3) 

and  the  fluid  force  normal  to  the  cable  axis  is: 

Eno(l  -h  So)  =  —pConDoU  COS  po  | Z7  COS  p^  1  v/TT^ 

where  Do  is  the  diameter  of  the  unstretched  cable,  p 
is  the  density  of  water,  and  U  is  the  current  speed  as 
a  function  of  depth. 

The  frictional  and  normal  drag  coefflcients  are 
represented  by  Cot  and  Con  respectively.  Cot  is  not 
important  for  static  calculations  since  the  moorings 
are,  for  the  most  part,  nearly  vertical  and  the  current 
,  is  assumed  to  be  horizontal.  A  value  of  Cot  =  0.003 
-  is  most  often  used  based  on  calculations  of  local  skin 
friction  for  steady  turbulent  flow  along  a  long  cylin¬ 
der  (White,  1974).  The  normal  drag  coefficient  of  a 
cable  is  influenced  by  vortex-induced  vibrations  and 
can  range  from  a  value  of  Con  =  1.2  for  a  nonvibrat¬ 
ing,  smooth  circular  cylinder  to  values  greater  than 
Con  ^  2.0  for  cables  vibrating  with  amplitudes  on  the 
order  of  one  cable  diameter.  Kim  et  d  (1986)  mea¬ 
sured  drag  coefficients  of  1.4  to  1.6  for  a  long  vertical 
cables  that  were  lightly  tensioned  and  hanging  in  a 
shear  current.  Grosenbaugh  et  d,  (1991)  measured 


drag  coefficients  for  long  vertical  tow  cables  that  had 
similcLT  tension  to  that  of  a  mooring  cable  and  found 
drag  coefficients  between  1.6  atnd  2.0  in  shear  currents. 
For  the  mooring  calculations  presented  in  this  paper, 


we  use  a  value  of  Con  =  1-6. 

There  are  also  two  geometric 

relations  that 

give  the  position  of  the  cable  in 
coordinates.  These  are: 

cartesian  x,  x- 

=  (1  -re^)  cos  a>0 
as 

(5) 

—  =  (1  rieo)sinac 

(6) 

The  two  cable  equations  and  the  two  geomet¬ 
ric  relations  along  with  appropriate  boundary  condi¬ 
tions  can  be  integrated  using  Runge-Kutta  methods 
(Forsyxhe  et  aL,  1977)  to  give  efiective  static  tension 
Toj  cable  angle  Oo,  and  r,  c-conhg<jration.  The  inte¬ 
gration  begins  at  the  suriace  buoy  (s  =  L  where  L  is 
the  total  unstretched  length  of  the  composite  mooring 
line)  with  boundary  conditions: 

Toil)  =  F6,[cosfflo(I)]-' 

XoiL)  =  0  (7) 

ro(L)  =  0 

Fho  is  t:he  steauiy  fluid  drag  force  acting  on  the  buoy 
and  is  given  by: 

F'bo  =  -pCob--^bUb\Ub\  (8) 

where  Cob  is  the  drag  coefficient,  is  the  horizontal 
projected  area  of  the  buoy  and  Ub  is  the  current  speed 
at  the  surface  buoy. 

These  boundary  conditions  assume  that  ail 
the  forces  act  through  the  center  of  gravity  and  that 
the  chainge  in  the  draft  of  the  buoy  due  to  the  cable 
tension  is  small.  More  complicated  boundary  condi¬ 
tions  that  allow  for  the  calculation  of  the  equilibrium 
draft  and  pitch  can  be  used.  However  we  find  them  to 
be  unnecessary  because  of  the  large  waterplanfe  area 
of  the  surface  buoys  used  on  oceano^aphic  moorings 
and  the  large  weight  of  attached  instruments  which 
helps  keep  the  buoy  nearly  vertical.  The  final  bound¬ 
ary  condition  is  that  the  bottom  of  the  mooring  cable 
at  s  =  0  is  located  precisely  at  the  anchor.  This  can 
be  expressed  mathematically  as: 

ZoiO)  =  -d  (9) 


w'here  d  is  the  water  depth. 

The  static  solution  is  found  by  guessing  the  ca¬ 
ble  angle  at  the  top  of  the  mooring  line,  using  the  guess 
to  calculate  the  other  initiad  conditions  from  equation 
(7),  and  then  integrating  the  cable  equations  to  find 
To,  do,  Xo,  and  Zo  as  a  function  of  s.  If  aquation  (9)  is 
satisfied,  the  calculated  values  represent  the  true  so¬ 
lution.  If  the  end  of  the  cable  is  above  or  below  the 
anchor,  then  the  calculations  must  be  repeated  with 
a  modified  value  of  do(T).  The  Newton- Raphson  it¬ 
eration  process  (Forsythe  et  al.,  1977)  continues  until 
the  correct  value  of  (Po{L)  is  chosen,  and  the  bottom 
boundary  constraint  is  satisfied. 

The  attached  instruments  are  treated  as  cable 
elements  and  assigned  an  equivalent  weight  in  water 
per  unit  length,  a  diameter,  and  drag  coefficients  that 
give  the  proper  forces.  This  is  made  necessary  by  the 
fact  that  some  instrument  cages  are  more  than  2  me¬ 
ters  long.  The  combined  instrument  string  can  repre¬ 
sent  29  to  40  meters  of  the  total  length  of  the  mooring 
line.  .A.n  alternate  (and  equivalent)  method  (Goodman 
et  al.  1972)  is  to  treat  the  instruments  as  point  masses 
though  the  extra  length  of  the  instruments  still  needs 
to  be  taken  into  account. 

In  analyzing  the  statics  of  an  oceanographic 
mooring,  the  standard  procedure  is  to  examine  two 
different  current  cases.  The  first  case,  called  the  “sur¬ 
vival  current^’,  represents  an  extreme  current  profile. 
Using  this  in  the  static  calculations  gives  the  maximum 
static  load  at  each  location  of  the  mooring  line.  It 
also  provides  the  input  in  terms  of  the  equilibrium  ca¬ 
ble  configuration  and  tension  for  dynamic  calculations 
which  are  needed  to  predict  total  extreme  loads.  The 
second  case,  called  the  “design  current*'  corresponds  to 
n’pical  operating  conditions.  The  results  of  static  cal¬ 
culations  using  this  profile  are  used  to  verify  that  the 
angles  of  the  instruments  are  within  the  design  range 
(<  15°  from  the  vertical).  The  “design”  configuration 
and  tensions  are  also  used  as  input  for  the  numerical 
codes  to  predict  the  dynamic  tensions  that  are  used  in 
the  fatigue  amalysis. 

Dynamic  Calculations 

The  d^Tiamics  are  calculated  by  linearizing  the 
motion  about  the  static  configuration.  We  assume  that 
all  the  dependent  variables  (i.e.  the  tension,  the  cable 
angle,  the  longitudinal  and  transverse  motion,  zind  the 
strain)  axe  made  up  of  a  static  part  and  a  small  dy¬ 
namic  part.  .A.S  shown  in  reference  (Triantafyllou  et  a/., 
1987),  the  two-dimensional  d>mamic  equations  for  a 
cable  element  of  mass  per  unit  unstretched  length  rrio, 
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longitudinal  added  mass  per  unit  unstretched  length  cransverse  harmonic  motion.  The  value  of  coefficient 
and  transverse  added  mass  per  unit  unstreteched  depends  on  the  current  velocity  and  is  given  by: 

length  TTiga  can  then  be  written  as: 


,  .  d^p  dTi  ddr, 

(mo+TTip^)— =  — -To— +Fa  (10) 


(mo 


.S-q_d6o  doi  ,  dT^  . 

(11) 


ds 


ds 


^  -  doodsq  =  Ti/EA  (12) 

—  -r  do^dsp  =  (1  ri  eo)  (13) 

where  p  is  the  amplitude  of  longitudinal  motion,  q  is 
the  amplitude  of  transverse  motion.  T  is  the  d^mamic 
tension,  and  Q\  is  the  dynamic  angle.  The  values  of  To 
and  6o  are  known  from  the  static  solution.  The  toted 
tension  and  angle  are  then  given  by: 

T  =  To^T,  (14) 

ri  di  (15) 

The  added  mass  in  the  longitudinal  direction  is  zero 
for  the  mooring  cable  but  has  a  significant  value  for 
the  attached  instruments. 

The  fluid  forces  are  approximated  by  equiva¬ 
lent  linear  terms  such  that: 


(17) 

We  determine  the  values  of  bp  and  by  assuming 
harmonic  motion  and  equating  energy  lost  per  cy¬ 
cle  for  the  linear  and  nonlinear  drag  representations 
(Paulling,  1979).  The  expressions  used  in  the  calcu¬ 
lations  which  take  into  ax:count  the  magnitude  of  the 
steady  current  axe  (Triantafyllou  et  al^  1987): 


bp  =  CtpCoi^du  IpI  (18) 

67  =  c^pConduj  1^1  (19) 


2 

=  - 
T 


(2^Dl) 


ri  UjiarcsinUR 


(20) 


where  Ur  is  the  absolute  value  of  the  ratio  berween  the 
current  velociu"  and  the  cable  velocity  either  \U/Ljp\  or 
\U/uq\.  For  values  of  Ur  >  1,  Ce  =  Ur, 

For  a  narrow  banded  random  input  such  as 
that  for  irregular  seas,  the  expression  bp  and  b^  become 
(Caughey,  1963): 


=  CtpCot^dffp  (21) 

bq  =  Dnd(Tq  (22) 

where  eXp  is  the  standard  deviation  of  the  longitudinal 
cable  veiociir/  and  cTp  the  standard  deviation  of  the 
transverse  cable  velocity.  The  coefficient  becomes: 

<=«  =  y  f  ^  (23) 

where  now  Ur  =  \JJ/cTp\  or  Ur  =  iCZ/a^l  depending 
on  the  direction  of  motion.  For  both  simple  harmonic 
motion  and  irregular  motion  defined  by  am  input  spec¬ 
trum,  the  damping  constants  bp  and  bq  depend  on  the 
solution  which  must  be  obtained  by  iteration. 

The  solution  in  the  frequency  domain  is  found 
by  assuming: 


Ti  =  (24) 

4>i  =  (25) 

P  =  (26) 

<7  =  (27) 

where  Ti,  (pi,  p,  and  q  are  complex  amplitudes.  The 

expressions  are  substituted  into  equations  (10-13)  to 
obtain  ordinary/  differ ential  equations  with  respect  to 
the  spatial  \-ariable  s.  We  can  write  these  in  matrix 
form  as: 


where  y  is  defined  as: 


f  fi\ 

P 

\  Q  J 


(28) 


(29) 


where  u  is  the  frequency  of  motion  and  |p|  and  |gj 
are  the  respective  amplitudes  of  the  longitudinal  and 
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M  is  the  mass  of  the  buoy,  /in,  /iss)  2.nd  /in  are 
frequency  dependent  added  mass  coefncients,  An,  A33, 
A55,  and  An  are  frequency  dependent  wave  damping 
coefficients,  C33  and  C55  are  the  restoring  constants 
for  heave  and  pitch,  and  Fr^i,  Fr,3,  Fr^s  a.re  the 
wave  forces  and  moments.  In  the  next  section,  we  de¬ 
scribe  how  to  calculate  the  hydrod^mamic  coefncients 
and  w'ave  forces  for  an  axisymmetric  buoy  and  give 
values  for  three  different  buoy  shapes. 

The  equations  of  motion  for  the  buoy  are  writ¬ 
ten  about  its  center  of  gravity.  Here,  we  allow  the 
connection  point  between  the  cable  and  the  buoy  to 
be  located  at  any  height  above  or  below  the  center  of 
gravity.  This  distance  is  represented  in  equation  (3  ^ ) 
by  r^. 

If  we  assume  that  the  buoy  motions  aire  har¬ 
monic  along  with  the  tension  and  cable  angle  at  the 
top  of  the  mooring,  then  equations  (35-37)  simplify  to: 


where  H  is  a  complex  matrix  containing  the  buoy  mass 
plus  hydrod^mamic  coefficients,  is  a  complex  vector 
which  represents  the  wave  forces,  and  G  is  complex 
matrix  that  is  a  function  of  the  static  tension  and  angle 
at  the  top  of  the  cable. 

The  boxmdar^'  condition  at  the  top  of  the  cable 
is  that  its  motion  is  compatible  with  the  motion  of  the 
buoy.  They  axe  related  by: 


cos  Oo  sin  <i>o 
—  sin  <po  cos  <po 


(39) 


where  R  is  a  2  x  3  matrix  that  relates  the  cable  co¬ 
ordinate  system  to  the  buoy  coordinated  system.  If 
we  left  multiply  both  sides  of  equation  (38)  by 
and  use  equation  (34)  to  eliminate  T^v^and  we  fi¬ 
nally  arrive  at  an  expression  for  the  motion  of  the  top 
of  the  mooring  cable  in  terms  of  known  input.  The 
expression  is: 


Once  p.v  a,nd  q,\r  are  found,  the  dynamic  tension  and 
dynamic  cable  angle  at  the  anchor  can  be  found  from: 


and  then  the  motion,  d\mamic  tension,  tind  d^mamic 
angle  at  nth  node  of  the  cable  can  calculated  from  the 
partial  matrix  product: 

(42) 

The  buoy  motion  is  found  by  inverting  equation  (39). 

For  harmonic  calculations  (i.e.  at  a  particu¬ 
lar  wave  frequency  and  wave  amplitude),  we  guess  chat 
the  amplitude  of  motion  at  each  node  pi  and  q{  is  equal 
to  the  wave  amplitude.  We  use  this  guess  to  determine 
the  B^’s  and  the  B{/s.  We  then  solve  equation  (40) 
for  Piv  and  qy  and  equation  (41)  for  Ti  and  oi.  Equa¬ 
tion  (42)  allows  us  to  solve  for  the  motion  along  the 
length  of  the  mooring  line  which  are  then  compared 
to  the  original  guess.  This  process  continues  until  the 
guess  and  solution  converge  within  a  user  specified  er¬ 
ror.  .4 way  from  the  natural  elastic  frequencies  of  the 
mooring  cable,  the  solution  process  converges  quickly, 
usually  in  two  iteration  steps.  Convergence  is  slower 
near  natural  frequencies  of  the  system. 

For  spectral  analysis  where  the  input  is  de¬ 
fined  by  a  wave  spectrum  St7(u),  the  standstrd  devi¬ 
ation  of  the  longitudinal  and  transverse  cable  veloci¬ 
ties  cTp  and  are  initiaLily  set  equal  to  the  standard 
deviation  of  the  wave  velocin".  This  is  used  to  calcu¬ 
late  the  transfer  functions  and  for 

the  longitudinal  and  transverse  mooring  motions.  The 
spectrom  of  the  cable  motions  Sp[uyS)  and  Sg(u,s) 
axe  then  found  from: 


Spiu.s)  =  \Hpiuj,s)\'S,iu)  (43) 

5,(:^,5)  =  (44) 

New  values  of  (Tp  and  are  found  by  integrating  the 
velocity"  spectra  which  axe  just  cj-  times  the  motion 
spectra.  The  calculations  are  then  repeated  until  con¬ 
vergence  is  achieved.  The  spectrum  of  the  tensions  at 
any  location  can  then  be  calculated  from: 


1  =  (40) 

IN  J  '■ 


Sr(y,5)  =lffT(w,5)|-S.,(L;)  (45) 
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Discus  Buoy 

Nominal  DiameTer  3.Q  m 

Waterplane  Area  6.2 

Displaced  Volume  2.1 

Metacenrric  Height  1. 1  m 
Moment  of  Inertia  1900 


Spherical  Buoy 

Nominal  Diameter  3.0  m 

Waterplane  Area  5.7 

Displaced  Volume  3.8  m^ 

Metacentric  Height  0.37  m 

Moment  of  Inertia  5200  kg-m^ 


Toroid  Buoy 

Nominal  Diameter  2.4  m 

Waterplane  Area  4.0  m^ 

Displaced  Volume  1.2  m^ 

Metacentric  Height  1.1  m 

Moment  of  Inertia  1400  kg-m^ 


Figure  S;  TT-picai  oceanographic  surface  buoys. 


Hydrodynamic  Coefficients  of  Oceanographic 
Buoys 

The  r^v'o-dimensional  numerical  method  just 
presented  incorporates  the  exciting  wave  forces,  mo¬ 
ments,  and  hydrod^mamic  coefficients  of  the  surface 
buoy  into  the  caicuiations.  Here  we  present  the  wave 
forces  and  moments  and  hydrodynamic  coefficients  for 
three  oceanographic  buoys  that  axe  currently  being 
used  on  surface  moorings  (Fig.  8).  The  S^foot  toroid 
buoy  was  used  on  the  original  oceanographic  surface 
moorings  of  the  1960*s,  and  is  still  used  today  on  AT¬ 
LAS  moorings  (Milbum  and  McClain,  1986).  The  3- 
meter  discus  buoy  was  developed  in  the  lace  1970*s  to 
fill  the  need  of  increased  payload  and  reserve  buoy¬ 
ancy  (Berteaux  and  Roy,  1986).  It  has  become*  the 
workhorse  of  the  oceanographic  community.  The  10- 
foot  spherical  buoy  was  developed  recently  as  a  buoy 
with  large  reserve  buoyancy  and  a  shape  that  has  im¬ 
proved  seakeeping  qualities.  To  date,  it  has  been  used 
exclusively  on  oceanographic  engineering  test  moor¬ 
ings.  Buoy  characteristics  are  listed  in  Table  1. 

All  of  the  buoys  in  Fig.  8  are  axisymmecric 


about  the  vertical  axis.  This  faot  is  used  to  simplify 
the  calculations.  We  consider  a  free  floating  vertical 
axisymmetric  body  excited  by  regular  waves  of  ampli¬ 
tude  Tf  and  frequenc^r  u  propagating  in  water  depth 
d.  Cylindrical  coordinates  (r,5,z)  are  introduced  with 
origin  at  the  undisturbed  water  surface  and  the  2-axis 
positive  upwards.  Viscous  effects  are  neglected  and 
the  assumption  is  made  that  the  fluid  is  incompress¬ 
ible.  The  motions  of  the  body  and  the  fluid  are  as¬ 
sumed  to  be  small  so  that  linear  potential  theory  can 
be  applied  to  describe  the  fluid  motion.  The  flow  is 
~  governed  by  the  velocity  potentiad  ^  which  satisfies 
"  the  Laplace  equation  within  the  fluid  domain  and  the 
appropriate  boundary’'  conditions  on  the  free  surface, 
the  sea  bottom,  the  wetted  surface  of  the  body,  and 
the  radiation  condition  at  infinity  (Mei,  1983). 

We  examine  the  in-plane  motion  (surge  fi, 
heave  <^3,  aind  pitch  ^5)  of  the  buoy  in  regular  waves. 
The  velocity  potential  for  this  case  can  be  expressed 
as: 

^ir,9,z,t)  =Rc{0{r,9,z)e-'“'^]  (45) 
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with 


unit  normal  vector  fi  and  position  vector  f  such  that: 


o{r,0,  z)  =  6o{r,9,  z)  ^  <h{r,d,=) -r  ^  ^j6{r,9,z) 

j= 1,3.3 

(46) 

Here,  is  the  velocity  potentiaii  of  the  undisturbed 
incident  harmonic  wave,  07  is  the  difiraction  potential 
for  the  body  fixed  in  the  wave,  and  o,-  is  the  radiation 
potential  resulting  from  the  forced  body  motion  in  the 
jth  mode  of  motion  with  unit  velodrr  amplitude 

The  veiocit}’'  potential  of  the  ^disturbed  inci¬ 
dent  wave  can  be  expressed  in  r/iindrical  coordinates 
as: 


6q{t,9,z)- 

.  cosh\k(  z^d)] 


/csinn  Ka 


where  Jm  is  the  mth  order  Bessel  function  of  first  kind 
and  €m  is  the  Neumann's  symbol  such  that  £0  =  1  and 
£rn  =  -  for  m  >  1.  The  frequency  c,'  and  wave  number 
k  are  related  by  the  dispersion  relation: 


cj-  =  gkzanhkd 


In  a  similar  manner,  the  diffraction  and  radiation  po¬ 
tentials  can  be  written  in  the  form: 


00 

(pj  (r,  9,  z)  =  -iuTi  ^  z)  cos  m9  (49) 

m=sO 

for  j  =  1,3, 5, 7.  In  this  case,  the  unknown  complex 
functions  ^jm  are  established  through  the  method  of 
matched  axisymmetric  eigenfunction  expansions.  This 
method  was  introduced  in  (Garrett,  1971)  for  the  solu¬ 
tion  of  the  diffraction  problem  around  a  vertical  cylin¬ 
der  and  used  for  both  the  diffraction  and  radiation 
problem  of  a  vertical  cylinder  in  (Yeung,  1981).  The 
technique  that  we  use  for  calculating  the  wave  exdting 
forces  and  hydro dmamic  coefn dents  is  an  extension  to 
arbitrarily  shaped  vertical  axisymmetric  bodies  that  is 
described  in  (Mavrakos,  1985;  Kokkmowrachos  et  ai, 
1986;  Mavrakos,  1988). 

Given  the  solutions  for  the  velodty  potentials 
<P7  and  dj,  we  can  determine  the  exdtifig  wave  forces 
and  moments  used  in  equations  (35-37)  from: 


Fr^j  =  -iup  J  y*  (do  ^  dr)  rijdS  (50) 

where  the  integration  is  taken  over  the  wetted  surface 
area  of  the  buoy.  Here,  Uj  is  denned  in  terms  of  the 


(nx,n2,n3)=n  (n^, ns , n^)  =  r  x  h  (51) 

where  we  are  interested  in  the  planer  components  ni, 
n3,  and  ns.  The  hydrodynamic  added  mass  and 
damping  Xij  are  given  by: 


jlij  -r  ^X{j  —  p 


o^UidS 


The  method  outlined  above  has  been  used  to 
obtain  the  hydrodymamic  chciract eristics  of  the  three 
oceanographic  buoys  depicted  in  Fig.  8.  The  excit¬ 
ing  forces,  moments,  and  hydrodymamic  coefficients  for 
surge,  heave,  and  pitch  are  given  in  Figs.  9-11.  The 
curves  correspond  to  d  =  100  m,  though  the  results 
are  nearly  identical  to  those  found  for  d  =  oo. 

.\s  can  be  seen  from  Fig.  9,  where  the  toroidal 
oceanographic  buoy  is  considered,  a  peculiar  behav¬ 
ior  of  both  the  heave  exciting  forces  and  the  hydrody¬ 
namic  parameters  is  obtained  near  /  =  0.7  Hz.  This 
behavior,  which  is  tymical  of  toroidal  bodies  (Garrett, 
1970;  Newman,  1977;  Miioh,  1983),  daises  from  the 
resonance  heave  oscillation  of  the  free  surface  enclosed 
by  the  toroid.  In  the  neighborhood  of  this  resonance, 
negative  values  of  the  added  mass  coefficients  are  ob¬ 
tained  though  the  damping  coefficients  always  remain 
positive. 


Sample  Calculations  of  Surface  Moorings 


Sample  caicuiacions  are  presented  to  show 
the  dymamic  characteristics  of  oceanographic  surface 
moorings.  We  will  use  the  design  for  the  1991  MLML 
mooring.  This  mooring  was  an  inverse-catenary  moor¬ 
ing  deployed  in  2850  m  of  water.  In  order  to  examine 
characteristics  of  taut  moorings,  we  perform  some  of 
the  calculations  using  a  depth  of  3700  m  which  is  equal 
to  the  total  length  of  the  mooring  line  (Fig.  12). 

The  mooring  consisted  of  a  90  m  long  in¬ 
strument  string  followed  by  500  m  of  7/ 16-inch  and 
1710  m  of  3/8-inch  diameter  wdre  rope,  680  m  of  7/8- 
inch  diameter  nylon  rope,  690  m  of  1-1/8-inch  di¬ 
ameter  polypropylene  rope,  and  130  m  of  glass-balls 
mounted  on  chain  along  with  the  anchor  tackle.  The 
instrument  string  had  five  Multi- Variable  Moored  Sys¬ 
tems  (MVMS)  which  are  modified  VMCMs  with  ad¬ 
ditional  sensors  mounted  in  the  cage,  four  bio-optical 
sensors  (BOMS)  which,  like  a  SEACAT,  these  con¬ 
sist  of  an  instrument  pressure  housing  that  is  mounted 
on  a  steel  rod,  one  .A.DCP,  one  engineering  instru¬ 
ment  that  is  attached  in-line  to  measure  mooring-line 


11 


Frequency  (Hz)  Frequency  (Hz) 

Figure  9:  Hydrodynamic  coefficients,  wave  forces,  and  wave  moments  for  toroid  buoy. 
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Figure  10:  Hydrodynamic  coefficients,  wave  forces,  and  wave  moments  for  discus  buoy. 
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Figure  11;  Hydrod^nrxamic  coemdents,  wave  forces,  and  wave  momenrs  for  spherical  buoy. 


tension,  and  one  cusrom-made  bioluminescence  sensor 
(MOORDEX).  Importani  characteristics  of  the  instru¬ 
ments  are  given  in  the  Appendix. 

The  static  configurations  are  calculated  as  de¬ 
scribed  at  the  beginning  of  this  section,  and  these  are 
used  as  input  for  the  dymamic  calculations.  We  use  a 
uniform  (^dth  depth)  current  of  0.02  m/s  for  slack  con¬ 
ditions.  For  moderate  currents  (i.e.  design  current), 
we  use  a  profile  that  is  equal  to  0.60  m/s  from  the 
surface  down  to  a  depth  of  100  m,  decreasing  linearly 
to  0.25  m/s  at  1500  m  depth,  and  then  decreasing  to 
0.05  m/s  at  the  anchor.  The  biggest  difference  between 
the  taut  mooring  and  the  inverse-catenary  mooring 
is  the  static  tension  in  moderate  currents  (Fig.  12). 
The  larger  tensions  in  the  taut  mooring  (30,400  N  for 
the  taut  mooring  versus  17,700  N  for  the  inverse  cate¬ 
nary  in  moderate  currents)  increase"  the  chances  of  the 
mooring  components  breaking,  the  buoy  being  pulled 
underwater,  or  having  its  anchor  drag. 

For  the  initial  dynamic  calculation,  we  deter¬ 
mine  the  transfer  function  of  the  dymamic  tension  at 
the  top  of  the  taut  mooring  with  a  3-meter  discus  buoy 
in  Sea  States  3  and  7.  We  use  the  static  configuration 
corresponding  to  the  taut  mooring  in  moderate  cur¬ 


rents.  The  key  features  of  the  transfer  functions  (Fig. 
13)  aus  the  low  frequency  beha\dor  corresponding  to 
the  acceleration  of  the  mass  of  the  instruments  and 
wire  rope  (plus  the  drag  of  the  instraments  in  high 
sea  states),  the  fundamental  elastic  natural  frequency 
of  the  nyion  rope  near  0.3  Hz  and  the  heave  natural 
frequency  of  the  buoy  near  0.4  Hz.  These  natural  fre¬ 
quencies  dre  more  damped  in  the  higher  sea  states. 
The  d>mamic  coefScient  of  friction  for  the  wire  rope 
and  synthetic  rope  is  taken  as  (7/  =  0.003  (Grosen- 
baugh,  1995a). 

Additional  resonant  frequencies  above  0.5  Hz 
are  present  but  aue  relatively  unimportant  because  of 
the  limited  wave  energy  at  these  frequencies.  This  can 
be  seen  in  plots  of  the  motion  spectra  of  the  buoy  for 
surge,  heave,  and,  pitch  (Fig.  14)  where,  for  the  most 
part,  the  buoy  follows  waves.  The  pitch  calculations 
are  based  on  potential  flow,  and  the  real  response  most 
likely  has  more  damping  which  would  lower  the  reso¬ 
nant  response. 

We  recalculate  the  transfer  functions  of  the 
dynamic  tension  using  the  the  8-foot  toroid  buoy  and 
10- foot  spherical  buoy  and  compare  the  results  of  ail 
three  buoys  in  Fig.  15.  The  discus  and  spherical 
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Figure  13:  Transfer  function  of  the  tension  at  the  top 
of  a  taut  mooring  with  a  discus  buoy  in  Sea  State  3 
and  Sea  State  7. 


Figure  15:  Comparison  of  transfer  functions  of  the  ten¬ 
sion  at  the  top  of  a  taut  mooring  with  three  different 
buoy  shapes. 
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Figure  16:  Comparison  of  transfer  functions  of  the 
tension  at  the  top  of  a  taut  mooring  and  an  inverse- 
catenary  mooring  in  moderate  and  slack  current  con¬ 
ditions.  The  analytical  model  represented  by  the  solid 
curve  is  described  in  Section  3. 


Sea  State 

Discus 

Sphere 

Toroid  | 

3 

2220  N 

2150  N 

2350  N 

1 

4130  N 

3950  N 

4250  N 

Table  1:  Comparison  of  the  standard  deviations  of  the 
tensions  at  the  top  of  taut  mooring  for  different  buoy 
shapes. 

buoys  show  little  difference  except  for  slight  offset  of 
the  buoy  natural  frequency.  The  toroid  buoy  is  similar 
at  low  frequencies  but  shows  large  differences  near  0.75 
Hz  where  the  free-surface  resonance  within  the  center 
of  the  toroid  occurs.  The  tensions  on  the  mooring  line 
calculated  from  Equation  (45)  show  differences  in  Sea 
State  3  and  in  Sea  State  7  with  the  spherical  buoy 
giving  the  best  result  (Table  1).  ^  . 

The  next  calculation  compares  the  dynamic 
tensions  of  a  taut  mooring  and  an  inverse-catenary 
moorings  in  slack  and  moderate  current  conditions. 
Each  mooring  uses  a  3-meter  discus  buoy.  In  a  weak 
current,  the  loop  in  the  synthetic  rope  of  the  inverse- 
catenary  mooring  is  present.  In  a  moderate  current, 
the  loop  is  pulled  nearly  straight,  and  the  total  moor¬ 
ing  line  forms  a  catenary  with  positive  curvature. 


Figure  17:  Comparison  between  measured  tension 
power  spectra  for  the  1991  MLML  mooring  and  pre¬ 
dictions  based  on  nximericai  simihation. 

The  d\mamic-tension  transfer  function  of  the  moor¬ 
ings  show  that  the  response  of  the  taut  and  inverse- 
catenary  moorings  in  moderate  currents  are  similar 
(Fig.  16).  In  slack  current  conditions  the  effects  of 
the  spring  in  the  synthetic  rope  are  absent  so  that  the 
resonance  at  0.3  Hz  is  not  present  and  the  heave  nat¬ 
ural  frequency  of  the  buoy  is  reduced.  The  solid  curve 
corresponds  to  an  analytical  model  that  is  described 
in  the  next  section. 

For  the  final  calculation  in  this  section,  we  use 
the  inverse-catenary  design  of  the  1991  MLML  moor¬ 
ing  and  compare  estimated  tension  spectra  for  six  dif- 


Sea  State 

Measured 

Simulation  | 

AnaJyticai  Model 

3 

I960  N 

2160  N 

2210  N 

4 

2340  N 

2520  N 

2570  N 

5 

2940  N 

3010  N 

3060  N 

5-H 

3530  N 

3690  N 

3810  N 

6 

3750  N 

3720  N 

3780  N 

7 

4020  N 

4050  N 

4120  N 

Table  2:  Comparison  of  the  standard  deviations  of  the 
tensions  at  the  top  of  1991  MLML  moorings  for  given 
sea  states. 
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ferenc  sea  states  with  actual  measurements  (Fig.  17). 
Qualitatively,  che  fit  to  the  data  is  good.  The  stan¬ 
dard  deviation  of  the  tensions,  given  in  Table  2,  show 
errors  of  11%  for  the  lowest  sea  state  and  errors  that 
are  less  than  3%  for  the  other  sea  states. 

ANALYTICAL  MODEL 

In  this  section,  we  present  an  analytical  model 
that  can  used  for  preliminary  d^mamic  analysis  of 
instrumented  oceanographic  surface  moorings.  The 
model  simplifies  the  problem  by  treating  only  the  ver¬ 
tical  motion  of  the  buoy  and  the  longitudinal  motion 
of  the  mooring  line  and  attached  instruments.  We  will 
show  through  comparison  with  the  two-dimensional 
numerical  model  that  the  anadytical  model  captures 
all  the  important  dynamic  effects  and  gives  accurate 
predictions  of  the  dynamic  tension  at  the  top  of  the 
mooring  line. 

The  analytical  model  vras  proposed  by 
Grosenbaugh  (1995b)  for  predicting  the  tensions  at  the 
top  of  a  mooring  line  where  the  tensions  are  great¬ 
est.  The  model  simplifies  the  problem  by  trearing 
only  the  vertical  motion  of  the  buoy  and  the  longi¬ 
tudinal  motion  of  the  mooring  line  and  attached  in- 
stroments.  The  justification  for  this  assumption  is 
based  on  the  effects  of  hydrodynamic  damping.  The 
transverse  motions  that  are  excited  at  the  top  by  the 
buoy’s  surge  and  pitch  motions  dissipate  quickly  as 
the3/  travel  down  the  mooring  line  and  have  little  effect 
on  the  dynamic  tension.  Conversely,  the  longitudinal 
motions  are  maintained  by  the  hydrod^mamic  damping 
all  the  way  down  the  mooring  line  even  for  the  inverse- 
catenar;/  mooring  where  the  slope  of  the  mooring  line 
can  become  negative  in  slack  current  conditions. 

The  major  simplifying  assumption  of  the 
model  is  that  the  dynamics  of  the  instruments  and 
their  connecting  hardware  can  be  combined  into  a  sin¬ 
gle  point  mass  and  damper  that  hangs  just  below  the 
buoy.  The  instruments  are  ni^icaily  located  within  the 
top  few  hundred  meters  of  the  mooring.  Our  assump¬ 
tion  is  that  the  whole  instrument  string  moves  as  a 
single  unit  with  the  same  amplitude  and  phase  of  the 
surface  buoy.  This  creates  an  effect  where  the  d3mamic 
tension  at  top  of  the  mooring  line  is  a 'function  of  the 
sum  of  their  mass,  added  mass,  and  viscous  damping. 
In  fact,  for  a  heavily  instrumented  mooring,  90-95%  of 
the  d}Tiamic  tension  can  be  attributed  to  the  instru¬ 
ments.  The  remaining  tension  comes  from  the  elastic 
strain  of  the  wire  and  synthetic  rope. 

Using  these  assumptions,  the  system  is  rep¬ 
resented  by  three  coupled  differencial  equations:  one 


describing  the  heave  motion  of  the  buoy,  one  describ¬ 
ing  the  vertical  motion  of  the  instruments  (i.e.  point 
mass  and  damper),  and  one  describing  the  longitu¬ 
dinal  motion  of  the  combination  wire  rope/synthetic 
rope  mooring  line.  These  can  be  solved  to  give  the 
transfer  function  of  the  tension  at  the  top  of  the  moor¬ 
ing.  The  spectral  response  for  a  given  sea  state  can 
then  be  calculated  from  equation  (45).  For  a  string 
of  instruments  writh  total  mass  and  added  mass  of 
including  any  hardware  used  to  connect  the  instru¬ 
ments  together  and  an  equivalent  linearized  damping 
constant  3  which  is  related  to  the  number  and  size  of 
the  instruments,  the  transfer  function  of  the  tension 
is: 


=  Hz{u)  [-.Mw-  -f  iScj-  (54) 

cot  tan/c^Lti/)] 

The  first  two  terms  in  brackets  account  for  in¬ 
ertia  and  drag  of  the  instrument  string  located  above 
the  wrire  rope.  The  drag  force  on  the  instruments  is 
a  function  of  the  standard  deviation  of  the  buoy^s 
heave  velocity  w'nich  is  calculated  from  the  buoy's 
heave  transfer  function  as  described  below.  The  third 
term  in  brcickets  is  associated  with  the  strain  of  the 
synthetic  rope  where  E,  is  Young's  modulus,  A,  is 
the  cross-sectional  area,  and  L,  is  the  length  of  the 
synthetic  rope.  This  term  becomes  important  when 
ksLs  =  t/2  which  defines  the  fundamental  elastic  nat¬ 
ural  frequencTr-.  The  wave  number  of  the  synthetic  rope 
ks  is  defined  by: 


T7T,a;'  — 

) 


(55) 


where  is  the  standard  deviation  of  the  mooring 
line's  longitudinal  velocity  averaged  over  the  length 
of  the  synthetic  rope,  m,  is  the  mass  per  unit  length, 
and  bs  is  the  equivalent  linearized  damping  constant 
..per  unit  length  of  the  synthetic  rope.  The  last  term 
in  brackets  is  associated  writh  the  response  of  the  wrire 
rope  where  is  Young's  modulus,  A^,  is  the  cross- 
sectional  area,  and  is  the  length  of  the  wrire  rope. 
The  wave  number  k^  is  given  by: 


k 


w 


a  .  , 

TTT^CJ*  —  KTpO^jjiJ 
Eiif  44u/ 


(56) 


where  m,^  is  the  mass  per  unit  length  and  b^,  is  the 
equivalent  linearized  damping  constant  per  unit  length 
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of  the  wire  rope.  .4gain,  o-p  is  averaged  over  the  length 
of  the  wire  rope. 

Finally,  the  term  Hz{^)  is  the  transfer  func¬ 
tion  of  the  buoy’s  heave  response  in  waves  taking  into 
account  the  presence  of  the  instruments  and  the  moor¬ 
ing  line.  It  is  given  by: 

Hz{^)  =  [C'ss  -r  (5/) 

cci^Es -A-ghg  CQZ  kgLis  “•  E^A.itfkxu  tan  | 

where  Mz  -  M  ^  Mzz  ^  M  ^  Ma  is  the  sum  of  the 
mass  and  auided  mass  of  the  buoy  and  the  instrument 
string,  Bz  ^  Bzz  ^  is  the  sum  of  the  damping  of  uhe 
buoy  and  the  instrument  string,  and 

a=  fl>  (58) 

The  heave  motion  of  the  buoy  for  a  given  sea  state  is 
given  bv: 

Sz{u)  =  iHt;(w)l-5„(w)  (59) 

where  Sz{'-^)  is  che  power  spectrum  of  the  buoy  mo- 
tion.  Integration  of  times  the  motion  spectrum 
gives  standaird  de'^dation  of  the  heave  velocity  which  is 
used  to  calculate  the  damping  coefficient  of  the  instru¬ 
ments. 

The  analytical  model  is  compared  in  Fig.  16 
to  the  three  different  configurations  of  the  1991  MLML 
inverse  catenary  mooring  shown  in  Fig.  12.  The  com¬ 
parison  between  the  analytical  model  and  the  numer¬ 
ical  model  is  made  using  the  transfer  function  of  the 
dynamic  tension  given  by  Equation  (55 j  and 

the  parameters  listed  in  Table  3. 

For  both  small  and  large  sea  states,  the  ana¬ 
lytical  model  gives  similar  restiits  to  those  of  the  nu¬ 
merical  model  of  the  taut  mooring  and  the  inverse- 
catenary  mooring  in  moderate  currents.  For  the  1991 
MLML  mooring,  the  configuration  of  the  mooring  line 
could  be  represented  by  the  moderate-current  case 
over  95%  of  the  time.  The  spectra  in  Fig.  17  all  cor¬ 
respond  to  moderate  to  strong  currents,  and  the  ana¬ 
lytical  model  gives  excellent  fits  to  the  full-scale  mea¬ 
surements  (Grosenbaugh,  1995b).  The  standard  devi¬ 
ations  of  the  tension  for  the  six  different  sea  states^  are 
nearly  identical  to  the  results  of  the  numerical  model 
(Table  2).  This  shows  the  validit}^  of  the  assumption 
that  the  heave  motion  of  the  buoy  and  the  longitudi¬ 
nal  motion  of  the  mooring  line  are  the  cause  of  the 
dynamic  tension. 

An  aLd vantage  of  the  analjmcai  model  is  that 
it  shows  explicitly  the  individual  effects  of  the  instru¬ 
ments,  mooring  hardware,  wire  rope  and  synthetic 


Instrument 

String 

Total  mass 

Total  added  mass 

Damping  constant  (5) 
Length 

1080  kg 

400  kg 

1260  .V  s-/m^ 

90  m 

Wire 

Rope 

Mass  per  length 

Young's  modulus  (E^) 
Cross-sectional  area  (A,*,) 
Length  (Lm;) 

Friction  coefficient  (C/) 

0.45 

1.0  X  10^^  N/m- 
1.2  X  10-*^  m2 
2110  m 

0.003 

Synthetic 

Rope 

Mass  per  length  (ms) 
Young's  modulus  {Eg) 
Cross-sectional  area  {Ag) 
Length  (Ls) 

Friction  Coemcient  (C/) 

0.33  kgm^'^ 

4,8  X  10^ 

5.0  X  10“^^  m2 
1370  m 

0.003 

Table  3;  Parameters  used  in  analytical  model  for  1991 
MLML  mooring. 

rope  that  go  into  the  overadl  dynamic  response.  This 
is  important  to  engineers  in  the  early  phases  of  the 
design  process  when  it  is  being  decided  how  many  in¬ 
struments  and  how  much  wire  and  synthetic  rope  will 
be  used  in  the  mooring.  The  anaiyticad  model  is  an 
order  of  magnitude  faster  than  numerical  simulations 
and  gives  accurate  results  for  moorings  with  the  instru¬ 
ments  concentrated  in  the  top  300  m  of  the  mooring 
(Grosenbaugh,  1995b).  At  a  minimum,  the  analytical 
model  can  be  used  to  verify  the  results  of  numerical 
simulations. 

ULTIMATE  AND  FATIGUE  ANALYSIS 

The  flow  chart  shown  in  Fig.  5  serves  as  a 
guide  in  performing  ultimate  and  fatigue  failure  anal- 
-  ysis  of  a  mooring.  The  first  step  in  both  processes 
*  involves  characterizing  the  environmental  forcing  in 
terms  of  the  expected  currents  and  waves.  These  are 
used  to  calculate  static  and  dynamic  tensions  which 
are  then  incorporated  into  the  failure  analysis.  The 
static  and  dmamic  results  are  then  combined  during 
the  ultimate  failure  analysis.  The  fatigue  failure  anal¬ 
ysis  is  depends  primarily  on  the  dmamic  results. 

Design  aind  survival  current  profiles  are  deter¬ 
mined  from  historical  oceanographic  data.  Many  times 
the  moorings  are  placed  in  locations  that  have  never 
been  instrumented  and  so  estimates  must  be  made. 
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Sea 

State 

No. 

Wind 

Speed 

(m 

Significant 
Wave  Height 

(m) 

Freq. 

(Kz) 

Percentage 
Probability 
of  Sea  State 

L 

0-3.5 

0.72 

0.204 

20.7 

2 

3. 5-5.5 

1.52 

0.149 

11.7 

3 

5. 5-8. 5 

1.93 

0.135 

10.5 

4 

8.5-11 

2.37 

0.125 

12.8 

5 

11-14 

3.05 

0.112 

17.5 

6 

14-17 

4.20 

0.099 

16.4 

17-20 

5.82 

0.087 

7.3 

S 

20-24 

7.20 

0.080 

2.2 

9 

>24 

9.00 

0.074 

0.4 

Table  4:  Sea  state  occurrences  for  Arabian  Sea  moor¬ 
ing  site. 


The  current  profiles  used  in  the  design  of  the  .Arabian 
Sea  surrace  mooring  were  determined  from  ship  obser¬ 
vations  of  surface  currents  and  estimates  of  how  they 
changed  with  depth. 

A  description  of  the  wave  environment  is  con¬ 
structed  from  observations  within  or  near  the  deploy¬ 
ment  area.  .A  source  of  such  data  is  the  U.S.  Navy 
Marine  Climatic  Atlas  of  the  World  (Meserve,  1974) 
which  can  be  used  to  construct  a  listing  of  the  signif¬ 
icant  wave  height,  peak  wave  frequencn'’,  and  proba- 
bilitv'  of  occurrence  for  a  range  of  sea  states.  Wave 
data  for  the  Arabian  Sea  surface  mooring  site  are 
presented  in  Table  4.  A  wave  amplitude  spectrum 
corresponding  to  each  sea  state  can  be  constructed 
by  assuming  that  the  shape  is  given  by  a  modified 
Pierson-Moskowitz  spectram  for  niily-deveioped  seas 
or  a  JONSW'AP  spectrum  for  fetch  limited  seas  (Fal- 
tisen,  1990). 

The  next  step  is  to  calculate  the  static  tension 
and  equilibrium  mooring  configuration  corresponding 
to  the  design  and  survi\*al  current  profiles.  The  con¬ 
figuration  is  then  used  to  determine  the  tension  power 
spectrjm  for  each  sea  state  of  interest.  This  can  be 
done  by  calculating  the  transfer  function  of  the  ten¬ 
sion  numerically  and  using  Equation  (45)  to  deter¬ 
mine  St{^)-  Alternatively,  the  analytical  model  *can 
be  used,  if  applicable,  to  calculate  Ht{o)  at  the  top  of 
the  mooring  where  the  tensions  are  the  greatest,  and 
this  can  be  used  for  the  components  of  the  instrument 
string. 

The  standard  deviation  of  the  tension  ar  is 
needed  for  both  the  ultimate  and  fatigue  failure  anal¬ 
ysis.  This  is  determined  from  the  tension  spectra  by 


X  10* 


Tima  (3| 


Figure  18:  Time  series  of  the  tension  at  the  top  of  the 
1991  MLML  mooring  in  Sea  State  7. 

computing: 

=  (^J  5T(u;)dj^  (60) 

In  addition  for  fatigue  failure  analysis,  the  mean  fre¬ 
quency/  (in  Hertz)  of  the  tension  fm  5or  each  sea  state 
is  needed  and  is  determined  from: 


Ultimate  Failure  Analysis 

Ultimate  failure  analysis  involves  determining 
the  maximum  expected  load  on  a  mooring  component 
and  comparing  this  value  to  the  ultimate  strength.  If 
the  extreme  load  is  greater  than  the  ultimate  strength 
then  the  part  is  likely  to  fail.  In  practice  for  oceano¬ 
graphic  moorings,  a  safety  factor  of  2.5  is  applied  so 
that  the  design  criterion  for  all  mooring  components 
is: 

Tr  ^  ^ultimate  //?n\ 

extreme  <  - -  (0.6) 

ib.O 

The  instantaneous  tension  on  any  mooring 
component  is  the  combination  of  the  static  tension  To 
and  the  dynamic  tension  Ti  as  given  by  Equation  (14). 
For  oceanographic  moorings  that  are  heaveiy  instru¬ 
mented,  the  dvmamic  tension  can  be  of  the  same  mag¬ 
nitude  as  the  static  tension  (Fig.  18)  and,  thus,  cannot 
be  ignored  in  calculating  the  total  extreme  load. 

The  static  component  of  the  extreme  tension 
77  is  found  using  the  survival  current  profile.  The 
contribution  from  the  dynamic  tension  is  detemined 
by  first  calculating  crx  using  the  survival  cxirrent  pro¬ 
file  for  the  equilibrium  mooring-line  configuration  and 
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the  amplitude  spectrum  of  the  largest  expected  sea 
state.  The  extreme  d\mamic  tension  is  then  based  on 
choosing  the  average  tension  of  the  highest  1/n  ten¬ 
sion  amplitudes.  We  write  this  value  as  where 

we  choose  n  =  1000.  For  the  largest  sea  state  in  Table 
5  (Sea  State  7)  which  occurs  0.1%  of  the  time  over  an 
eight  month  depio\nnent,  turns  out  to  be  the 

average  of  the  five  highest  dynamic  tensions  assuming 
a  mean  wave  nequency  of  0.25  Hz. 

Water  wave  amplitudes  closely  follow  a 
Rayleigh  distribution.  If  we  assume  that  the  ampli¬ 
tudes  of  the  tension  also  follow  a  Rayleigh  distribution, 
we  can  write  the  extreme  d}mamic  tension  as: 


where  F  is  a  gamma  function  and  (Tt  is  the  standard 
deviation  of  the  tension  corresponding  to  a  particu¬ 
lar  sea  state.  The  parameters  q  and  y  are  fatigue 
constants  related  to  the  material  and  geometry  of  the 
given  mooring  component. 

The  material  constants  are  found  from  labora¬ 
tory  experiments  in  which  a  number  of  identical  moor¬ 
ing  components  are  c^'cied  with  constant-amplitude 
dynamic  tension  until  they  fail.  The  number  of  cycles 
at  which  the  failure  occurs  is  assumed  to  be  related  to 
the  amplitude  of  the  d>mamic  tension  through: 

(66) 


2.1/1000  _  (63) 

where  again  (jf  is  calculated  from  the  surh\ni  cur¬ 
rent  and  the  largest  expected  sea  state.  The  extreme 
tension  is  then: 


r«rrcm.  (64) 

This  formula  assumes  the  sundvai  current  occurs  at 
the  same  time  as  the  maximum  d^mamic  tension,  and 
thus  provides  additional  margin  of  safec%^  in  the  calcu¬ 
lation. 

Fatigue  Failure  Analysis 

To  predict  the  fatigue  life  of  a  mooring  com¬ 
ponent,  we  adopt  the  Palmgren- Miner  theory  for  cu¬ 
mulative  damage  (Miner,  1945).  The  theory  is  linear 
and  assumes  that  fatigue  damage  occurring  from  a  sin¬ 
gle  CTr'cle  of  loading  at  one  amplitude  can  be  added  to 
the  sum  of  the  damage  from  ail  previous  cycles  even  if 
the  previous  amplitudes  were  different.  If  a  mooring 
component  will  break  after  a  toted  of  N  cycles  when 
subjected  to  a  constant  amplitude  d^mamic  tension  T, 
then  it  follows  from  the  Palmgren-Miner  theory  that 
the  fraction  of  damage  after  one  cycle  is 

If  the  same  part  is  subjected  to  a  narrow- 
banded  random  loading  such  as  that  caused  by  wave- 
induced  motions  during  a  particular  sea  state,  then^the 
fraction  of  damage  over  one  cycle  is  the  expected  v^ue 
of  iV“^.  If  we  assTime  again  that  the  peak  tensions  fol¬ 
low  a  Rayleigh  probability  distribution,  then  it  is  easy 
to  show  (Triantafvdlou  et  ai.,  1986)  that  the  expected 
value  of  corresponding  to  a  given  sea  state  is 

£:(,v-t)  =  r(i4.i)  (65) 


Equation  (66)  is  used  widely  to  approximace  compo¬ 
nents  made  of  metal  including  wire  rope.  Synthetic 
rope  can  aiso  be  approximated  with  this  equation  but 
with  a  lesser  degree  of  accuracr/.  By  performing  tests 
with  different  ^•aiues  of  T,  we  can  determine  q  and  y 
using  least-squares  techniques. 

We  determined  material  constants  for  various 
sizes  of  galvanized-steel  sling  links  and  shackles  like 
the  ones  used  on  the  MLML  moorings  by  choosing  12 
parts  cind  cv'ciing  them  until  one  part  failed.  By  doing 
this  at  different  loads,  we  found  that  the  links  can  be 
represented  with  g  =  3.3  and  y  =  l.Ty^  (Fig.  19)  and 
the  shackles  can  be  represented  by  g  =  3.7  and  y  = 
l-^Vu  (Fig.  20)  where  is  the  ultimate  strength  of  the 
given  component.  This  method  of  taking  the  extreme 
minimum  -v-alue  was  chosen  because  our  emphasis  is  on 
insuring  the  survivability  of  the  mooring. 

We  aiso  tested  the  wire  rope  with  swage  fit¬ 
tings.  This  data  is  presented  in  Fig.  21  along  with 
other  wire-rope  fatigue  data  from.  Here  the  material 
constants  are  g  =  4.6  and  y  =  2.4yu. 

The  material  values  for  the  links  and  shackles 
correspond  to  tests  performed  at  a  mean  tension  of 
17,800  N.  To  account  for  differences  in  the  me3xi  load, 
we  determine  an  effective  y  and  use  this  in  Equation 
^,(65).  The  effective  value  y’  is  given  by  (Collins,  1993): 


yu  —  17,  800  \ 

Vu  ""  'Pm  ) 


(67) 


where  Tm  is  the  mean  tension  which  is  determined 
from  the  static  analyis. 

The  amount  of  fatigue  damage  due  to  a  given 
sea  state  is  the  fraction  of  damage  chat  occurs  over 
one  cycle  times  the  total  number  of  cycles  chat  occur 
during  chat  sea  state.  For  a  deployment  that  lasts  a 
total  of  r  seconds,  the  number  of  cycles  u  chat  occur 
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for  a  given  sea  3T:ate  is: 


Figure  19:  Cyclic  fatigue  curve  for  links. 


Figure  20:  Cyclic  fatigue  curve  for  shackles. 


^  /p)  (68) 

where  /p)  is  the  percentage  probability  of  occur¬ 
rence  for  a  given  sea  state  with  significant  wave  height 
/i,  and  peak  frequency  /p  (Table  4).  The  fatig'je  dam¬ 
age  for  a  given  sea  state  is  chen: 


D^uE{N^^)  (69) 


The  total  amount  of  fatigue  damage  for  the  whole  de¬ 
ployment  is  the  sum  of  the  damage  over  all  sea  states. 
Thus,  the  criteria  for  predicting  a  fatigue  mlure  is: 


_ 1 _ 

factor  of  safety 


(70) 


where  Dj  is  the  damage  associated  with  the  jth  sea 
state. 

The  factor  of  safety  accounts  for  the  uncer¬ 
tainties  caused  by  the  order  of  the  loading  and  by 
the  effects  of  corrosion.  For  quasi-random  loading,  the 
Miner's  sum  given  by  the  left  hand  side  of  Equation 
(70)  ranges  between  0.6  and  1.6.  Corrosion,  on  the 
other  hand,  can  reduce  fatigue  strength  by  as  much  as 
half  (Collins,  1993).  Accordingly,  we  use  a  factor  of 
safety  of  4.0,  accepting  a  mooring  component  only  if 
the  analysis  shows  that  it  loses  no  more  than  a  quarter 
of  its  life  during  a  deployment. 

This  analysis  was  applied  in  retrospect  to 
the  1989  and  1991  MLML  moorings  (Grosenbaugh, 
1995c).  In  particular,  it  predicted  that  the  5/8-inch 
pear-shaped  sling  link  shown  in  Fig.  4  would  last  44.4 
days  compared  to  the  actual  time  to  failure  of  70  days. 
In  addition,  the  5/ 8-inch  shackles  used  at  the  top  of 
the  mooring  line  would  have  had  a  fractional  loss-of-life 
of  1.19.  Based  on  our  criteria  of  safety  factor  of  four, 
this  part  would  not  have  been  acceptable.  The  fatigue 
resistance  of  the  1991  mooring  was  improved  by  using 
an  inverse  catenary  design  which  helped  to  lower  the 
^mean  tension,  using  7/8-inch  sling  links  and  3/4-inch 
“shackles  in  the  upper  part  of  the  mooring,  and  reduc¬ 
ing  the  deployment  time  by  15%.  The  7/S-inch  sling 
links  were  predicted  to  have  a  fractional  loss-of-life  of 
0.19  while  the  shackles  were  predicted  to  have  a  frac¬ 
tional  loss  of  0.23.  The  1991  mooring  was  recovered 
with  the  mooring  hardwaire  showing  little  wear. 


Design  Example 


Figure  21;  Cyclic  fatigue  curve  for  wire  rope. 
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We  will  analyse  two  components  of  the  1995 
Arabian  Sea  Surface  Mooring  to  demonstrate  the  de¬ 
sign  techniques  for  determining  ultimate  strength  and 


fatigue  resistance.  The  components  are  a  3 /4-inch 
chain  shackle  that  is  located  at  :he  top  of  the  moor¬ 
ing  where  tensions  are  the  highest  and  the  3/S-inch 
wire  rope  which  begins  at  the  bottom  of  the  instru¬ 
ment  string.  The  steps  in  this  exaimpie  are  the  same 
for  ail  other  mooring  components  so  this  can  ser^^e  a 
step-by-step  guide  for  engineering  analysis. 

The  Arabian  Sea  surface  mooring  was  an 
inverse- cat enar}'  mooring  (Fig.  22)  containing  approx¬ 
imately  1150  m  of  1  1/S-inch  diameter  polypropylene 
rope,  1680  m  of  7/8-mch  diameter  nylon  rope,  and 
1610  m  of  3/8-inch  and  7/ 16-inch  diameter  wire  rope. 
.Above  this  was  a  250  m  long  instrument  string  con¬ 
taining  1/2-inch  diameter  wire  rope,  3/4-inch  chain, 
and  a  suite  of  instruments  that  included  5  V2vICMs, 
four  MVMSs,  four  SEACATs,  one  large  custom  made 
instrument  for  detecting  zooplankton  with  sound  {Hol¬ 
iday  InstTument) ,  and  13  temperature  sensors. 

The  purpose  of  the  mooring  was  to  study 
ocean  surface  mixing  in  a  region  where  cooling  may 
occur  during  the  summer.  In  the  .Arabian  Sea,  it 
is  thought  that  strong  monsoon  winds  during  June, 
July,  and  August  may  cause  increased  mixing  with  the 
cooler  w’ater  from  below  the  surface.  This  hypothesis 
could  only  be  answered  with  a  surface  mooring. 

As  seen  in  Table  4,  the  Arabian  Sea  is  char¬ 
acterized  by  strong  winds  and  high  sea  states.  During 
10%  of  the  deployment,  the  winds  are  expected  to  be 
gale  force.  During  the  monsoon  months  this  percent¬ 
age  increases  to  25%.  With  the  large  number  of  instru¬ 
ments  on  the  mooring,  d^mamic  tensions  will  be  high. 
Currents  are  also  strong  creating  large  static  tensions. 

The  3/4-inch  shackles  are  used  to  attach  ail  of 
the  instruments  to  the  mooring  line,  the  most  heaveiy 
tensioned  shackle  is  located  just  above  the  top  VMCM. 
Thus,  for  tensions  in  the  auialysis,  we  will  use  the  ten¬ 
sion  at  the  top  of  the  mooring  line  for  the  calculations. 
Fig.  23  shows  the  mooring  configuration  and  static 
tension  at  the  buoy  in  the  design  and  survi\'al  cur¬ 
rents.  These  were  calculated  using  the  static  analysis 
procedures  outlined  in  Section  2, 

The  d^mamic  tensions  are  calculated  for  each 
of  the  sea  states  in  Table  4.  We  wdll  use  the  numerical 
model  for  a  random  sea  input,  relying  on  equations 
(43-45)  ,to  calculate  the  tension  power  spectra  of  the 
final  mooring  design  shown  in  Fig.  22.  For  the  wave 
spectra,  we  use  modified  Pierson-Moskowitz  spectra. 
This  is  conservative  as  fetch-limited  spectra  may  be 
more  appropriate.  The  tension  spectra  calculated  in 
this  manner  are  given  in  Fig.  24  and  the  standard 
de'viations  of  the  tension  and  average  frequency  are 


Figure  22:  Schematic  drawing  of  the  1995  Arabian  Sea 
mooring. 
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Figure  23;  Stacie  design  and  survival  configurations 
for  the  Arabian  Sea  mooring. 


Figure  24;  Tension  power  spectra  at  the  top  of  the 
Arabiain  Sea  Mooring  corresponding  to  the  sea  states 
listed  in  Table  5. 


Sea  State 

Design  Current 
<7T  (N)  /m  (Hz) 

Survival  Current 
cjf  (N) 

1 

2240 

0.391 

— 

2 

2880 

0.359 

3 

3100 

0.349 

— 

4 

3340 

0.342 

— 

5 

3570 

0.332 

— 

6 

3970 

0.321 

— 

7 

4410 

0.308 

— 

8 

4760 

0.299 

— 

9 

5250 

0.288 

5110 

Table  5;  Standard  deviation  and  mean  frequency  of 
tensions  at  the  top  of  the  1995  Arabian  Sea  surface 
mooring. 

given  in  Table  5.  We  note  that  much  of  the  preiim- 
inaxy  d^mamic  analysis  of  the  design  of  the  Arabian 
Sea  mooring  was  performed  with  the  analytical  model 
as  the  assumptions  of  nearly  taut  mooring  line  and 
instrument'String  length  less  than  300  m  apply. 

The  extreme  dynamic  tension  is  calculated  us¬ 
ing  Equation  (63)  with  the  extreme  standard  devia¬ 
tion  of  o-f  =  5110  N  corresponding  to  Sea  State  9 
and  the  su]rri\-al  current  profile  (last  column  of  Ta¬ 
ble  5).  We  then  determine  the  total  extreme  tension 
using  Equation  (64)  with  the  extreme  static  tension 
of  27  =  28, 100  N  which  corresponds  to  the  survival 
current  (Fig.  23).  Thus, 

=  23, 100  W  ^  3.S5  X  5110 //  (71) 

=  52,700  iV 

The  ultimate  strength  of  a  3/4-inch  shackle  is  254,000 
N.  This  represents  a  safety  factor  of  5.3  which  easily 
satisfies  the  design  constraints. 

The  fatigue  analysis  covers  a  six-month  de¬ 
ployment  (184  days).  This  corresponds  to  16  million 
seconds  or  roughly  five  million  cycles  of  loading.  Using 
the  probabilities  in  Tables  4,  the  average  frequencies 
in  Table  5  (corresponding  to  the  design  current),  and 
Equation  (68),  we  calculate  the  number  of  en'des  corre¬ 
sponding  to  each  sea  state  (column  1  in  Table  6).  The 
standard  deviations  in  Table  5  along  with  the  mate¬ 
rial  fatigue  constants  represented  by  the  dashed  line  in 
Fig.  20  (y=1.5  cuud  q=3.7)  are  substituted  into  Equa¬ 
tion  (65)  to  give  the  damage  per  cycle  for  each  sea 
state.  The  fraction  of  life  lost  due  to  the  jth  sea  state 
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Sea  State 

Cycles 

Galvanizec 
Fraction  Lost 

Shot-Peened 
Fraction  Lost 

l.2sr  x!.J“ 

0.043 

O.OLO 

2 

0.668  xlO® 

0.059 

0.013 

3 

0.583  xi0“ 

0.067 

0.015 

4 

0.696  xlO® 

0.106 

0.024 

5 

0.924  xlO®  ' 

o 

OO 

o 

0.040 

6 

0.337  xl0° 

0.242 

0.054 

7 

0.382  xlO^ 

0.163 

0.036 

S 

0.105  xl0° 

0.059 

0.013 

g 

0.013  xl0° 

0.015 

0,003 

1 

t 

1  Sum  =  0.936 

1  Sum  =  0.208 

Table  6:  Summary  of  fatigue  analysis  for  3/4-inch  gal¬ 
vanized  and  shot-peened  shackles. 

is  given  by  Equation  (69)  (column  3  in  Table  6).  Sum¬ 
ming  :hese  values  gives  zhe  total  fraction  of  life  lost 
during  the  deployment.  The  value  for  the  gai'vanized 
shackle  is  0.936  which  gives  a  safety  factor  of  only  1.1 
which  is  below  the  acceptable  value  of  4.0. 

Going  to  a  larger  shackle  was  not  an  option  be¬ 
cause  the  shackle  pin  would  not  ht  through  the  ends 
of  the  instruments  and  it  would  be  too  expensive  to 
modir/  all  of  the  instrument  cages.  The  solution  in 
this  case  was  to  use  a  shackle  that  w^as  shot  peened  and 
painted  with  an  corrosion- resistant  coating.  The  ulti¬ 
mate  strength  of  this  part  is  the  same  as  the  galvanized 
component.  The  fatigue  resistance  of  the  shot-peened 
shackle,  estimated  from  the  data  point  in  Fig.  20,  is 
increased  by  50%.  Using  fatigue  constants  of  (y=2.25 
and  q=3.7),  we  recalculated  the  damage  per  cycle  and 
the  fraction  of  life  lost  due  to  the  jth  sea  state  (column 
4  in  Table  6).  Summing  these  numbers  gives  a  'value 
for  the  shot-peened  shackle  of  0.208  which  produces 
an  acceptable  safety  factor  of  4.8. 

To  aLnai>'ze  the  strength  and  fatigue  properties 
of  the  3/8-inch  ware  rope,  we  take  the  static  cuid  dy¬ 
namic  tensions  at  the  junction  between  the  wire  rope 
and  the  instrument  string.  In  the  numerical  simulation 
this  involves  outputting  the  tension  at  the  intersect¬ 
ing  node.  The  extreme  static  tension  corresponding  to 
the  sur.’iv’al  current  is  20,300  N.  The  extreme  dynamic 
tension  corresponding  to  Sea  State  9  (and  the  surviv^- 
current  mooring  configuration)  is  found  from  the  nu¬ 
merically  calculated  \alue  for  the  standard  deviation 
of  the  tension  (last  column  of  Table  7).  The  extreme 
total  tension  is  then  found  from  Equation  (64): 

T,ztremc  =  20, 300  iV  ^  3.85  X  1420  iV  (72) 


Sea  State 

Design  Current 
ar  (N)  (Hz) 

Survival  Current 
erf  (N) 

1 

720 

0.504 

— 

2 

360 

0.466 

— 

3 

910 

0.454 

— 

4 

970 

0.445 

— 

5 

1020 

0.432 

— 

6 

1120 

0.417 

— 

7 

1230 

0.400 

— 

3 

1310 

0.388 

— 

9 

1430 

0.374 

1420 

Table  7;  Standard  deviation  and  mean  frequency  of 
tensions  in  3/8-Inch  Wire  Rope  1995  Arabian  Sea  Sur¬ 
face  Mooring. 


Sea  Scate 

Cycles 

Fraction  Lost 

1 

1.659  xlO® 

0.0004 

2 

0.867  xlO® 

0.0004 

3 

0.753  xl0« 

0.0005 

4 

0.906  xlO® 

0.0008 

5 

1.202  xl0« 

0.0013 

6 

1.087  xlO® 

0.0019 

1 

0.496  xlO® 

0.0013 

8 

0.136  xlO® 

0.0005 

9 

0.024  xlO® 

0.0001 

Sum  Fracxion  Lost  =  0.0072 


Table  8:  Summary  of  fatigue  analysis  for  3 /S- inch  wire 
rope. 

=  25,800  iV 

In  ultimate  strength  tests  of  3/8-Lnch  wire  rope  with 
swage  fittings  on  the  end,  the  wire  breaks  a  few  cen¬ 
timeters  away  from  the  swage  joint  with  an  average 
strength  of  65,900  N.  This  represents  a  safety  factor 
of  2.6  which  satisfies  the  design  constraints.  Three 
hundred  meters  of  3/8-inch  wire  rope  were  eliminated 
from  the  original  mooring  design  to  achieve  a  safety 
factor  above  2.5. 

For  fatigue  analysis,  we  use  the  probabilities 
in  Tables  4,  the  average  frequencies  in  Table  7  (corre¬ 
sponding  to  the  design  current),  and  Equation  (68)  to 
determine  the  number  of  cycles  corresponding  to  each 
sea  state  (column  1  in  Table  8).  The  standard  devia¬ 
tions  in  Table  7  (for  the  design  current)  along  with  the 
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material  fatigue  constants  represented  by  the  dashed 
line  in  Fig.  21  (y=2.4  and  q=4.6)  are  substituted  into 
equation  (65)  to  give  the  damage  per  crrcie  for  each  sea 
state.  The  fraction  of  life  lost  due  to  the  ;th  sea  state 
is  given  by  equation  (69),  Summing  these  values  gives 
a  total  fraction  of  life  lost  during  the  deployment  of 
only  0.0072  (Table  8).  Thus,  fatigue  of  the  wire  rope 
is  not  a  concern. 

The  analysis  of  the  shackle  and  the  wire  rope 
are  examples  where  ultimate  strength  analysis  (of 
the  wire  rope)  and  fatigue  analysis  (of  the  shackle) 
led  to  design  changes.  Similar  analyses  were  per¬ 
formed  on  all  mooring  components  including  the  in¬ 
strument  strength  members.  Fatigue  results  led  to 
using  oval-shaped  links  (instead  of  pear-shaped  links) 
and  strengthening  the  welds  at  the  ends  of  the  VMCM 
and  MVMS  instrument  cages. 


CONCLUSIONS 

Dynamic  analysis  is  key  to  designing  oceano¬ 
graphic  surface  moorings  for  rough  water  emdron- 
ments.  It  is  important  in  predicting  both  extreme 
loads  and  fatigue  damage.  We  have  presented  a  fre¬ 
quency  domain  numerical  method  and  an  analytical 
model  that  allow  designers  to  perform  these  calcula¬ 
tions  in  different  sea  states  and  have  shown  how  to 
incorporate  the  results  into  ultimate  strength  and  fa¬ 
tigue  analysis. 

The  anah-tical  model  described  in  Section  3  is 
useful  in  the  early  stages  of  design  because  it  shows  ex¬ 
plicitly  the  individual  effects  of  the  instruments,  moor¬ 
ing  hardware,  wire  rope  and  synthetic  rope  that  go  into 
the  overall  dynamic  response.  It  is  orders  of  magnitude 
faster  than  numerical  simulations,  and  gives  accurate 
results  for  moorings  with  the  instruments  concentrated 
in  the  top  portion  of  the  mooring. 

Hydrod}mamic  coefficients  and  wave  exciting 
forces  of  different  shaped  ocemographic  buoys  axe 
given  along  with  the  material  and  fatigue  parime- 
ters  for  the  components  and  instruments  that  make 
up  the  mooring  line  and  instrument  strings.  These  are 
used  with  the  results  from  static  and  d>TicLmic  analy¬ 
sis  to  predict  if  a  mooring  will  surrive  a  deployment. 
The  case  study  of  the  Arabian  Sea  surface  mooring 
shows  how  the  design  techniques  presented  in  this  pa¬ 
per  can  isolate  weaknesses  and  lead  to  improve  designs 
in  terms  of  ultimate  strength  and  fatigue  resistance. 
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APPENDIX 

The  foiloudng  tables  give  important  parame¬ 
ters  of  oceanographic  instruments  used  on  the  1991 
Marine- Light  Mixed-Layer  mooring  and  the  1995  Ara¬ 
bian  Sea  mooring.  These  values  are  used  for  the  static 
and  dynamic  calculations  presented  in  this  paper. 

In  the  tables.  Type  corresponds  to  the  the 
method  of  attachment.  Open  Cage  uses  a  strucrural 
frame  of  3 /4-inch  stadniess-steel  rod  with  the  instru¬ 
ment  suspended  inside  the  frame.  Examples  of  this 
are  the  VMCM  and  ADCP  shown  in  Fig.  2.  Side- 
Mounted  on  Rod  corresponds  to  the  instrument  hous¬ 
ing  clamped  to  a  single  piece  of  1-1/4-inch  diameter 
stainless-steel  rod  which  acts  as  the  strength  mem¬ 
ber.  .An  example  of  this  is  the  SEACAT  shown  in  Fig. 
2.  In-Line  Cylinder  corresponds  to  a  cylinder  pres¬ 
sure  housing  with  attachment  points  mounted  on  the 
end  caps  so  chat  pressure  housing  becomes  the  struc¬ 
tural  member  of  the  mooring  line.  .An  example  of  this 
is  the  engineering  instrument  that  measures  mooring¬ 
line  tension.  Clamp  to  Mooring  Line  uses  a  damp 
that  hold  the  instrument  in  place  on  the  mooring  line 
by  friction. 

The  added  mass  eind  drag  coeffidents  are  es¬ 
timates  made  by  the  authors.  The  drag  coeffident  is 
assumed  to  be  the  same  for  both  vertical  and  trans¬ 
verse  motion. 

Table  A1 


Vector  Measuring  Current  Meter  (VMCM) 

Type 

Open  Cage 

Length 

2.9  m 

Mass 

81  kg 

Weight  in  'Water 

495  N 

Vertical  .Added  Mass 

20  kg 

Transverse  .Added  Mass 

33  kg 

Vertical  Projected  Area 

0.09  m- 

Transverse  Projected  Area 

0.56  wr 

Drag  Coefficient 

1.0 
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Table  A2 

Multi- Variable  Moored  Syste 

m  (MVMS) 

Type 

Open  Cage 

Length 

2.9  m 

Mass 

94.6  kg 

Weight  in  Water 

581  N 

Vertical  Added  Mass 

21  kg 

Transverse  Added  Mass 

35  kg 

Vertical  Projected  Area 

0.10  m" 

Transverse  Projected  Area 

0.68  m- 

Drag  Coefficient 

1.0 

Table  A3 

Acousric  Doppler  Current  Profiler  (ADCP) 

Type 

Open  Cage 

Length 

2.3  m 

Mass 

107  kg 

Weight  in  Water 

516  N 

Vertical  Added  Mass 

27  kg 

Transverse  Added  Mass 

54  kg 

Verticsd  Projected  Area 

0.26  m- 

Transverse  Projected  Area 

0.93  m- 

Drag  Coefficient 

1.0 

Table  A4 

SEACAT 

Type  Side-Mounted  on  Rod 

Length 

2.0  m 

Mass 

18  kg 

Weight  in  Water 

112  N 

Vertical  Added  Mass 

0.5  kg 

Transverse  Added  Mass 

6.6  kg 

Vertical  Projected  Area 

0.01  m- 

Transverse  Projected  Area 

0.05  m- 

Drag  Coefficient 

1.0 

Table  A5 

Bio-Opticai  Moored  System  (BOMS) 

Type  Side-Mounted  on  Rod 

Length 

2.9  m 

Mass 

61.5  kg  \ 

Weight  in  Water 

'305  N 

Vertical  Added  Mass 

4.5  kg 

Transverse  Added  Mass 

25  kg 

Vertical  Projected  Area 

0.03  m- 

Transverse  Projected  Area 

0.15  m- 

Drag  Coefficient 

1.0 

Table  A  6 

Custom-Made  Bio-Optical  Sensor 

(MOORDEX) 

Type  Custom-Made  Attachment 

Length 

2.0  m 

Mass 

157  kg 

Weight  in  Water 

979  N 

Vertical  Added  Mass 

51  kg 

Transverse  Added  Mass 

32  kg 

Vertical  Projected  Area 

0.29  m- 

Transverse  Projected  Area 

0.19  mr 

Drag  Coefficient 

0.8 

Table  A7 

Holliday  Bio- Acoustic  Instrument 

Type 

Open  Cage 

Length 

2.9  m 

Mass 

138  kg 

Weight  in  Water 

1030  .N 

Vertical  Added  Mass 

20  kg 

Transverse  Added  Mass 

33  kg 

Vertical  Projected  Area 

0.09  m- 

Transverse  Projected  Area 

0.56  m" 

Drag  Coefficient 

1.0 

Table  AS 

Temperature  Pod  Sensors  (TPod) 

Type  Clamp  to  Mooring  Line 

Length 

0.5  m 

Mass 

4.5  kg 

Weight  in  Water 

35.6  N 

Vertical  Added  Mass 

0.6  kg 

Transverse  Added  Mass 

4.0  kg 

Vertical  Projected  Area 

0.01  m* 

Transverse  Projected  Area 

0.05  m* 

Drag  Coefficient 

1.0 

Table  A9 

Engineering  Instrument 

Type 

In-Line  Cylinder 

Length 

1.0  m 

Mass 

40  kg 

Weight  in  W*ater 

225  N 

Verticad  Added  Mass 

2.0  kg 

Transverse  Added  Mass 

19  kg 

Vertical  Projected  Area 

0.02  wr 

Transverse  Projected  Area 

0.08  m- 

Drag  Coefficient 

1.0 
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